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Abstract 

A new dark energy model in anisotropic Bianchi type-Ill space-time 
with variable equation of state (EoS) parameter has been investigated in 
the present paper. To get the deterministic model, we consider that the 
expansion 8 in the model is proportional to the eigen value a\ of the 
shear tensor a 3 i . The EoS parameter ui is found to be time dependent 
and its existing range for this model is in good agreement with the recent 
observations of SNe la data (Knop et al. 2003) and SNe la data with 
CMBR anisotropy and galaxy clustering statistics (Tegmark et al. 2004). 
It has been suggested that the dark energy that explains the observed 
accelerating expansion of the universe may arise due to the contribution 
to the vacuum energy of the EoS in a time dependent background. Some 
physical aspects of dark energy model are also discussed. 

Keywords : Bianchi-III universe, Dark energy, Variable EoS parameter 
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1 Introduction 

Recent cosmological observations obtained by SNe la (Garnavich et al. 1998a, 
1998b; Perlmutter et al. 1997, 1998, 1999; Riess et al. 1998; Schmidt et al. 
1998; Knop et al. 2003), WMAP (Bennett et al. 2003; Spergel et al. 2003), 
SDSS (Tegmark et al. 2004; Seljak et al. 2005; Adelman-McCarthy et al. 2006), 
Chandra X-ray observatory (Allen et al. 2004) indicate that the observable uni- 
verse is undergoing an accelerating expansion. To explain the cosmic positive 
acceleration, mysterious dark energy has been proposed. There are several dark 
energy models which can be distinguished by, for instance, their variable equa- 
tion of state (EoS) oj(t) = ^ (p i s the fluid pressure and p its energy density) 
during the evolution of the universe. By now, methods allowing for restoration 
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of the quantity uj(t) from expressional data have been developed (Sahni and 
Starobinsky 2006), and an analysis of the experimental data has been conducted 
to determine this parameter as a function of cosmological time (see Sahni et al. 
2008 and references therein). DE has been conventionally characterized by EoS 
parameter mentioned above which is not necessarily constant. Recently, the 
parameter uo{t) is calculated with some reasoning which reduced to some simple 
parametrization of the dependences by some authors (Huterer and Turner 2001; 
Weller and Albrecht 2002; Chevallier and Polarski 2001; Krauss et al. 2007; 
Usmani et al. 2008; Chen et al. 2009). The simplest DE candidate is the vac- 
uum energy (uj = —1), which is mathematically equivalent to the cosmological 
constant A (Allen et al. 2004; Sahni and Starobinsky 2000; Sola and Stcfan- 
cic 2005; Shapiro and Sola 2009). The other conventional alternatives, which 
can be described by minimally coupled scalar fields, are quintessence (uj > — 1) 
(Ratra and Peebles 1988; Wcttcrich 1988; Liddle and Scherrer 1999), phantom 
energy (uj < —1) (Caldwell 2002; Caldwell et al. 2003) and the combination 
of quintessence and phantom in a unified model, namely quintom (Feng et al. 
2005; Guo et al. 2005) as evolved and have time dependent EoS parameter. 
Some other limits obtained from observational results coming from SNe la data 
(Knop et al. 2003) and SNe la data collaborated with CMBR anisotropy and 
galaxy clustering statistics (Tegmark et al. 2004) are —1.67 < uj < —0.62 and 
— 1.33 < uj < —0.79 respectively. The latest results, obtained after a com- 
bination of cosmological datasets coming from CMB anisotropies, luminosity 
distances of high redshift type la supernovae and galaxy clustering, constrain 
the dark energy EoS to —1.44 < uj < —0.92 at 68% confidence level (Hinshaw 
et al. 2009; Komatsu et al. 2009). However, it is not at all obligatory to use a 
constant value of uj. Due to lack of observational evidence in making a distinc- 
tion between constant and variable uj, usually the equation of state parameter 
is considered as a constant (Kujat et al. 2002, Bartelmann et al. 2005) with 
phase wise value — 1,0, — | and +1 for vacuum fluid, dust fluid, radiation and 
stiff dominated universe, respectively. But in general, uj is a function of time 
or redshift (Jimenez 2003; Das et al. 2005; Ratra and Peebles 1988). For in- 
stance, quintessence models involving scalar fields give rise to time dependent 
EoS parameter uj (Turner and White 1997; Caldwell et al. 1998; Liddle and 
Scherrer 1999; Steinhardt et al. 1999). Some literature are also available on 
models with varying fields, such as cosmological models with variable EoS pa- 
rameter in Kaluza-Klein metric and wormholes (Rahaman et al. 2006, 2009). 
In recent years various form of time dependent uj have been used for variable A 
models (Mukhopadhyay et al. 2008, 2009; Usmani et al. 2008). Recently Ray 
et al. (2010), Mukhopadhyay et al. (2010), Akarsu and Kilinc (2010), Yadav 
(2010), Yadav & Yadav (2010), Pradhan et al. (2010c) and Kumar (2010) have 
obtained dark energy models with variable EoS parameter in different contexts. 

Cosmologists have proposed many candidates for dark energy to fit the cur- 
rent observations such as cosmological constant, tachyon, quintessence, phantom 
and so on. The major difference among these models are that they predict dif- 
ferent equation of state of the dark energy and different history of the cosmos 
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expansion. The simplest dark energy (DE) candidate is the cosmological con- 
stant A, but it needs some fine-tuning to satisfy the current value of the DE. 
Overduin and Cooperstock (1998), Sahni and Starobinsky (2000), Komatsu et 
al. (2009) have suggested some dynamic models, where A varies slowly with cos- 
mic time (t). Srivastava (2005), Jackiw (2000), Bertolami et al. (2004), Bento 
et al. (2002); Bilic et al. (2002) and Avclino et al. (2003) have considered 
Chaplygin gas and generalized Chaplygin gas as possible dark energy sources 
due to negative pressure. Other than these approaches, some authors have con- 
sidered modified gravitational action by adding a function f(R) (R being the 
Ricci scalar curvature) to Einstcin-Hilbert Lagrangian, where f(R) provides a 
gravitational alternative for DE causing late-time acceleration of the universe 
(Capozziello 2002; Caroll et al. 2004; Dolgov and Kawasaki 2003; Nojiri and 
Odintsov 2003, 2004; Abdalaa et al. 2005; Mena et al. 2006). Recently Gupta 
and Pradhan (2010) have presented an entirely new approach as cosmological 
nuclear energy is a possible candidate for DE. For detail informations regard- 
ing the dynamics of dark energy, the readers are advised to see the reviews 
by Copeland et al. 2006 and Nojiri and Odintsov (2007). The aforementioned 
models offer satisfactory description of dark-energy bahaviour and its observ- 
able feathers. In spite of the success of these attempts, the nature of DE is one 
of the greatest challenge of modern cosmology. 

Spatially homogeneous and anisotropic cosmological models play a signifi- 
cant role in the description of large scale behaviour of universe and such models 
have been widely studied in framework of General Relativity in search of a real- 
istic picture of the universe in its early stages. Yadav et al. (2007), Pradhan et 
al. (2010a, 2010b) have recently studied homogeneous and anisotropic Bianchi 
type-Ill space-time in context of massive strings. Recently Yadav (2010) has 
obtained Bianchi type-Ill anisotropic DE models with constant deceleration pa- 
rameter. In this paper, we have investigated a new anisotropic Bianchi type-Ill 
DE model with variable ui without assuming constant deceleration parameter. 
The out line of the paper is as follows: In Section 2, the metric and the field 
equations are described. Section 3 deals with the solution of the field equations. 
In Section 4, some physical aspects of the derived DE model are given. Finally, 
conclusions are summarized in the last Section 5. 

2 The Metric and Field Equations 

We consider the space-time of general Bianchi-III type with the metric 

ds 2 = -dt 2 + A 2 (t)dx 2 + B 2 (t)e- 2ax dy 2 + C 2 {t)dz 2 , (1) 
where a is constant. 

The simplest generalization of EoS parameter of perfect fluid may be to de- 
termine the EoS parameter separately on each spatial axis by preserving the 
diagonal form of the energy momentum tensor in a consistent way with the 
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considered metric. Therefore, the energy momentum tensor of fluid is taken as 

Tf - diag[7%,T},T2,T$]. (2) 
Thus, one may parameterize it as follows, 

T( = diag[p, -p x , -p y , -p z ] = diag[l, -w x , -uj y , -u z ]p 

— diag[X, — u, — (cv + 8), — (w + j)]p. (3) 

Here p is the energy density of fluid; p x , p y and p z are the pressures and ui x , 
oj y and uj z are the directional EoS parameters along the x, y and z axes respec- 
tively, u is the deviation-free EoS parameter of the fluid. We have parameterized 
the deviation from isotropy by setting uj x — to and then introducing skewness 
parameter S and 7 that are the deviations from oj along y and z axis respectively. 

The Einstein's field equations (in gravitational units c = 1, 8ttG = 1) read 
as 

4-\M=n, (4) 

where R\ is the Ricci tensor; R — g 1 -* Rij is the Ricci scalar. In a co-moving 
co-ordinate system, the Einstein's field equation (|4]) with ([3]) for the metric (fT]) 
subsequently lead to the following system of equations: 

r r. uti 

-up, (6) 

u + S)p, (7) 
a 2 , , 

= 0. (9) 

Here and in what follows an over dot denotes ordinary differentiation with re- 
spect to t. 

The spatial volume for the model (QJ is given by 

V 3 = ABCe~ ax . (10) 

We define a = {ABCe~ ax )^ as the average scale factor so that the Hubblc's 
parameter is anisotropic models may be defined as 
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We define the generalized mean Hubble's parameter H as 



H=-{H 1 +H 2 + H 3 ) 7 



(12) 



where Hi = ^, H 2 — 5 and H3 = 5 are the directional Hubble's parameters 
in the directions of x, y and z respectively. 

An important observational quantity is the deceleration parameter q, which 
is defined as 

<? = - — • (13) 

a z 

The scalar expansion 9, components of shear cr^ and the average anisotropy 
parameter Am are defined by 
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where AH = H - H(i = 1, 2, 3). 



3 Solutions of the Field Equations 

The field equations JS])-© are a system of five equations with seven unknown 
parameters A, B, C, p, uj, S and 7. Two additional constraints relating these 
parameters are required to obtain explicit solutions of the system. We assume 
that the expansion (6) in the model is proportional to the eigen value a 2 2 of the 
shear tensor o J i . This condition leads to 

B = £ 1 (AC) m \ (21) 
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where l\ and toi are arbitrary constants. The motive behind assuming this 
condition is explained with reference to Thorne (1967), the observations of the 
velocity-red-shift relation for extragalactic sources suggest that Hubble expan- 
sion of the universe is isotropic today within ss 30 per cent (Kantowski and 
Sachs 1966; Kristian and Sachs 1966). To put more precisely, red-shift studies 
place the limit 

^<0.3, 

on the ratio of shear a to Hubble constant H in the neighbourhood of our Galaxy 
today. Collins et al. (1980) have pointed out that for spatially homogeneous 
metric, the normal congruence to the homogeneous expansion satisfies that the 
condition j is constant. 

Equations © leads to 

A = mB, (22) 

where to is a positive integrating constant. We have revisited the solution re- 
cently obtained by Pradhan et al. (2010a). 

Using the value of Eq. ([22]) in {7J and subtract the result from equation {6j , we 
obtain the skewness parameter on y-axis is null i.e. 5 = 0. In this case Eqs. ([6]) 
and {Jj are reduced to 



B A BC AC 
B ~ A + ~BC ~ AC 



0. (23) 



Using {22J in {23]), we obtain 



, B BC\ 
(1-m )( 1 j + W =0. (24) 



As m ^ 0, Eq. {2J| gives 



B BC 
~B + ~BC 



= 0, (25) 



which on integration reduces to 



BC = k u (26) 



where k\ is an integrating constant. 
From Eqs. {2T) and ([22]). we obtain 



B = i 2 C\ (27) 



where £ 2 = l[- mi m e , I = Using {27]) in {25]), we get 



<*>-£, (28) 
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which on integration gives 

C= (£+l)fTT 



w 2 t+k > 



(29) 



where k 2 is an integrating constant. Using (|2"9")l in (|27|) and (j2"2"j) we obtain 



and 



b = e 2 (e + 



A = ml 2 (£+ 



£T 2 t + k > 



— 



u 2 



t + ko 



(30) 



(31) 



respectively. 

Hence the metric ((T|) reduces to the form 



ds 2 = -dt 2 



m£ 2 (£+l)TTT [j^t + k 2 



t -i 2 

7+T 



dx 1 



dy 2 
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(32) 



Using the suitable transformation 

m£ 2 (£ + l)7TTx = X, 
£ 2 (£+l)ihy = Y, 

(£ + 1)7TTZ = Z, 

££,2 ' 

the metric (f3"2")l is reduced to 

rfs 2 = -p 2 dT 2 + T 2L dX 2 + T 2L e~^ x dY 2 + T^dZ 2 

where 

M = (<? + 1)^TT, 

AT = m£ 2 M 1 
I 



(33) 



(34) 



L = 



(35) 
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4 Some Physical and Geometric Properties of 
the Model 

The expressions for the scalar of expansion 6, magnitude of shear a 2 , the average 
anisotropy parameter A m , deceleration parameter q and proper volume V for 
the DE model are given by 

= ^±^, (36) 
\T\ 2 

(37) 
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3 V P-PT 

*=-(2£Tir (39) 

iV 2 1 LC21+X) 

V = — M7T- L t- L . (40) 

771 

The rate of expansion Hi in the direction of x, y and z are given by 

*i=tf 9 = ^, (41) 

^3 = ^. (42) 
Hence the average generalized Hubble's parameter is given by 

H= L ^ + 1 \ (43) 

Using equations (|29|) - (j3Tj) in JS| , the energy density of the fluid is obtained as 

1 lfi 2 T 2 N 2 T 2L ' y ' 

Now, by using equations (|29|) . (|30]) and (|44)) in ((6|), the equation of state pa- 
rameter uj is obtained as 

L(t 2 +1+1) -1(1+1) 
U W Li(lW)' { ' 

Using equations (|3ll|) . (j3"T|) . f4"H) and f4"5j) in (|5|), the skewness parameter 7 (i.e. 
deviation from uj along z-axis) is derived as 

j^r + v£wW - 1) - L{U 2 - * - 1)] 
? = a2 ■ ( 4fi ) 

£/3 2 T 2 N 2 T 2L 
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T 



Figure 1: The plot of EoS parameter uj versus T and L 

From equation (|45|) . it is observed that the equation of state parameter uj is 
time dependent, it can be function of redshift z or scale factor a as well. The 
redshift dependence of uj can be linear like 

uj(z) = uj q + lj'z, (47) 

with J = (ff ) 2 = (see Refs. Huterer and Turner 2001; Weller and Albrecht 
2002) or nonlinear as 

u(z) = u + -j—, 48 
1 + z 

(Polarski and Chavellier 2001; Linder 2003). So, as for as the scale factor 
dependence of uj is concern, the parametrization is given by 

w(a) = ujo + w tt (l - a), (49) 

where ujq is the present value (a = 1) and oj a is the measure of the time variation 
J (Linder 2008). 

The SNe la data suggests that -1.67 < u < -0.62 (Knop et al. 2003) while 
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'■° 10.0 

T 



Figure 2: The plot of energy density p versus T and L 



the limit imposed on w by a combination of SNe la data with CMB anisotropy 
and galaxy clustering statistics is —1.33 < cj < —0.79 (Tegmark et al. 2004). 
So, if the present work is compared with experimental results mentioned above, 
then one can conclude that the limit of uj provided by equation (|45[) may ac- 
commodated with the acceptable range of EoS parameter. Also it is observed 
that either for T = or for mi = 0, the to vanishes and our model represents a 
dusty universe. 

For the value of ui to be in consistent with observation (Knop et al. 2003), 
we have the following general condition 

Ti < T < T 2 , (50) 

where 

_ 0.79£ap 
1 ~~ [n^/L{£(£ + 1) - i(0.38£ 2 - 0.24J? + 1)} 



(51) 
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and 



T 9 



N^fL{£{£ + 1) - L(0.67£ 2 - 2.34£ + 1)} 



(52) 



For this constrain, we obtain —1.67 < lo < —0.62, which is in good agreement 
with the limit obtained from observational results coming from SNe la data 
(Knop et al. 2003). For a special case for which £ = 1, a = 0.5, f3 = 2, 
L = 0.5, N = 1, where 0.899700 < T < 1.153226, we obtain the same limit 
-1.67 < lo < -0.62. 



From Eq. (|45j) . we have observed that, at cosmic time 

T= f lap 

\_Ny/L{U + W+T) - 1} 

lo = — 1 (i.e. cosmological constant dominated universe) and when 

[ £aB 

T< H 

_Ny/L{L£ + £(£+l) - 1} 

lo > — 1 (i.e. quintessence) and when 

£aB 

T> H 

Ny/L{L£ + £{£ + 1) - 1} 



(53) 



(54) 



(55) 



lu < — 1 (i.e. super quintessence or phantom fluid dominated universe) (Cald- 
well 2002). 

The variation of EoS parameter lo with cosmic time T is clearly shown in Fig- 
ures 1, as a representative case with appropriate choice of constants of inte- 
gration and other physical parameters using reasonably well known situations. 
From Figure 1, we conclude that in early stage of evolution of the universe, the 
EoS parameter lo was very small but positive (i.e, the universe was matter dom- 
inated) and at late time it is evolving with negative value (i.e. at the present 
time). The earlier real matter later on converted to the dark energy dominated 
phase of the universe. 



From Eq. ([44|). we note that p(t) is a decreasing function of time and p > for 
all times. This behaviour is clearly depicted in Figures 2 as a representative 
case with appropriate choice of constants of integration and other physical pa- 
rameters using reasonably well known situations. 

In absence of any curvature, matter energy density (f2 m ) and dark energy (ft\) 
are related by the equation 

n m + n A = i, (56) 



li 



where Q m = and fl\ = ^nW. Thus, equation (|55|) reduces to 



p A 



3iP 37J 2 



1. (57) 



Using equations (j43|) and (j44|) . in equation (|57|) . the cosmological constant is 
obtained as 

£ 2 (2£+l) 2 _ £ 2 (l + 2) a 2 
3^ 2 T 2 ^ 2 T 2 + N 2 T 2L ' [ ' 

In recent time the A-term has interested theoreticians and observers for vari- 



Lambda 




Figure 3: The plot of cosmological term A versus T and L 

ous reasons. The nontrivial role of the vacuum in the early universe generate 
a A-term that leads to inflationary phase. Observationally, this term provides 
an additional parameter to accommodate conflicting data on the values of the 
Hubble constant, the deceleration parameter, the density parameter and the 
age of the universe (for example, see the references Gunn and Tinsley 1975; 
Wampler and Burke 1988). The behaviour of the universe in this model will 
be determined by the cosmological term A ; this term has the same effect as a 
uniform mass density p e f / = —A, which is constant in time. A positive value 
of A corresponds to a negative effective mass density (repulsion). Hence, we 
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expect that in the universe with a positive value of A, the expansion will tend 
to accelerate; whereas in the universe with negative value of A, the expansion 
will slow down, stop and reverse. From Eq. (1551) . we see that the cosmological 
term A is a decreasing function of time and it approaches a small positive value 
at late time. From Figure 3, we note this behaviour of cosmological term A 
in the model. Recent cosmological observations suggest the existence of a pos- 
itive cosmological constant A with the magnitude A(Gh/c 3 ) « 10~ 123 . These 
observations on magnitude and red-shift of type la supernova suggest that our 
universe may be an accelerating one with induced cosmological density through 
the cosmological A-term. It is remarkable to mention here that the dark energy 
that explains the observed accelerating expansion of the universe may arise due 
to the contribution to the vacuum energy of the EoS in a time dependent back- 
ground. Thus, our DE model is consistent with the results of recent observations. 

From the above results, it can be seen that the spatial volume is zero at T = 
and it increases with the increase of T. This shows that the universe starts evolv- 
ing with zero volume at T = and expands with cosmic time T. From equations 
(|4Tj) and (|42|) . we observe that all the three directional Hubble parameters are 
zero at T — > oo. In derived model, the energy density tend to infinity at T = 0. 
The model has the point-type singularity at T = (MacCallum 1971). The 
shear scalar diverses at T = 0. As T — > oo, the scale factors A(t), B(t) and C(t) 
tend to infinity. The energy density becomes zero as T — > oo. The expansion 
scalar and shear scalar all tend to zero as T — > oo. The mean anisotropy param- 
eter are uniform throughout whole expansion of the universe when I ^ — \ but 
for I = — i it tends to infinity. This shows that the universe is expanding with 
the increase of cosmic time but the rate of expansion and shear scalar decrease 
to zero and tend to isotropic. At the initial stage of expansion, when p is large, 
the Hubble parameter is also large and with the expansion of the universe H, 
9 decrease as does p. Since |y = constant provided I ^ — ± the model does 
not approach isotropy at any time. The cosmological evolution of Bianchi type- 
Ill space-time is expansionary, with all the three scale factors monotonically 
increasing function of time. The dynamics of the mean anisotropy parameter 
depends on the value of I. 

From (|39|) we observe that 

(*) for I < --, q>0 

i.e., the model is decelerating and 

(ii) for £ > — -, q < 

i.e., the model is accelerating. Thus this case implies an accelerating model of 
the universe. Recent observations of type la supernovae (see Perlmutter et al., 
1999; Riess et al., 1998 and references therein) reveal that the present universe 
is in accelerating phase and deceleration parameter lies somewhere in the range 
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— 1 < q < 0. It follows that our DE model of the universe is consistent with the 
recent observations. 

5 Concluding Remarks 

An anisotropic Bianchi type-Ill DE model with variable EoS parameter lj has 
been investigated which is new and different from the other author's solutions. 
In the derived model, u> is obtained as time varying which is consistent with re- 
cent observations (Knop et al. 2003; Tegmark et al. 2004). It is observed that, 
in early stage, the equation of state parameter lj is positive i.e. the universe 
was matter dominated in early stage but in late time, the universe is evolving 
with negative values i.e the present epoch (see, Figure 1). Our DE model is in 
accelerating phase which is consistent with the recent observations. Thus the 
model (|34p represents a realistic model. 

In the derived DE model of the universe, the cosmological term is a de- 
creasing function of time and it approaches a small positive value at late time 
(i.e., the present epoch). The values of cosmological "constant" for the model 
is found to be small and positive, which is supported by the recent observations 
(Garnavich et al. 1998a, 1998b; Perlmutter et al. 1997, 1998, 1999; Riess et al. 
1998, 2000, 2004; Schmidt et al. 1998). 

The DE model is based on exact solution of Einstein's field equations for the 
anisotropic Bianchi-III space-time filled with perfect fluid. To my knowledge, 
the literature has hardly witnessed this sort of exact solution for anisotropic 
Bianchi-III space-time. So the derived DE model adds one more feather to the 
literature. 

The DE model presents the dynamics of EoS parameter lj provided by Eq. 
(p4"5)) may accommodated with the acceptable range —1.67 < lj < —0.62 of SNe 
la data (Knop et al. 2003). It is already observed and shown in previous sec- 
tion that for different cosmic times, we obtain cosmological constant dominated 
universe, quintessence and phantom fluid dominated universe (Caldwell 2002), 
representing the different phases of the universe through out the evolving pro- 
cess. Therefore, we can not rule out the possibility of anisotropic nature of DE 
at least in the framework of Bianchi-III space-time. 

Though there are many suspects (candidates) such as cosmological constant, 
vacuum energy, scalar field, brane world, cosmological nuclear-energy, etc. as 
reported in the vast literature for DE, the proposed model in this paper at least 
presents a new candidate (EoS parameter) as a possible suspect for the DE. 
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